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ABSTRACT 

We present an Ha spectral observation of a Balmer-dominated shock on the eastern side of Tycho's 
supernova remnant using SUBARU Telescope. Utilizing the High Dispersion Spectrograph (HDS), we 
measure the spatial variation of the line profile between preshock and postshock gas. Our observation 
clearly shows a broadening and centroid shift of the narrow-component postshock Ha line relative to 
the Ha emission from the preshock gas. The observation supports the existence of a thin precursor 
where gas is heated and accelerated ahead of the shock. Furthermore, the spatial profile of the emission 
ahead of the Balmer filament shows a gradual gradient in the Ha intensity and line width ahead of 
the shock. We propose that this region (~ 10 16 cm) is likely to be the spatially resolved precursor. 
The line width increases from ~ 30 kms^ 1 up to ~ 45 kms -1 and its central velocity shows a redshift 
of ~ 5 kms -1 across the shock front. The characteristics of the precursor are consistent with a cosmic 
ray precursor, although a possibility of a fast neutral precursor is not ruled out. 

Subject headings: ISM: supernova remnants - ISM: individual (Tycho, G120.1+1.4) - Shock Waves - 
line: profiles 



1. INTRODUCTION 

Balmer-dominated filaments are the signature of 
non-radiative shocks propagating into pa rtially neu- 
tral medium (Chevalier & Ray mond! Il978f) . The Ha 
line profile is composed of two distinctive components 
(narrow and broad) representing the velocity distri- 
bution of preshock a nd postshock gases, respectively 
()Chevalier et al.ll980h . High resolution spectroscopic ob- 
servations of several of these shocks have revealed that 
the width of the narrow component is unusually large 
(30 ~ 50 kms -1 ) for ambient neutral hydrogen, and it 
was proposed that the gas was heated in a precursor 
thin enough (< 10 17 cm) to avoid complete ionization 
of hy drogen (see I Ghavamian et al.ll200ltlSollerman et al.l 
120031 and references therein). 6 Two likely c andidates are 
cosm i c-ray (CR) and fast neutral precurso rs ( Smit h etlrtl 
119941: iHester et al.|[l99l iGhavamian et al.ll2001[ ). Both 
scenarios predict significant Doppler shifts of preshock 
gas. No clear indication of such a shift of the Ha nar- 
row component is reported (but see lLee et al.ll200~4l ). A 
careful comparison of postshock and preshock line pro- 
files with high spectral resolution is crucial for confirming 
the existence of such a precursor. 

The blast wave of Tycho's SNR, the historical 
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remnant of the 1572 supernova, ha s been known 
to ex h ibit Balmer-dominated em i ssion (|van den Berghl 
I197H iKamper fe van den Berghl 119781 ). The re- 
gion of the brightest Ha em ission (Knot g from 
IKamper fe van den Berghl [l978) is located along the 
north-eastern edge of the remnant. Faint, diffuse optical 
emission e xtends ~1 pc ahead of the Balmer-dominated 
filaments (G havamian et aill2000L GOO hereafter). This 
feature has been identified as a photoionization precur- 
sor (PIP) produced by photoionization of the preshock 
gas by He II A304A emission from behind the blast 
wave (GOO). The estimated temperature of the PIP is 
~ 1.2 x 10 4 K, which is not high enough to explain the 
observed width of the Balmer narrow component. 

The existence of Ha emission from the PIP makes 
Tycho an unique target for the study of nature of the 
thin precursor, as we can investigate the change of the 
line profile across this precursor by comparing the line 
emission from the PIP (preshock) and that of Knot g 
(postshock). In this Letter, we present high resolution 
(Echelle) long slit Ha spectra of Tycho Knot g and its 
PIP using the SUBA RU High Dispersion Spectrograph 
(|Noguchi et al.ll2002f ). Our observations reveal the line 
broadening and the Doppler shift of the the narrow com- 
ponent, providing strong evidence for the existence of a 
cosmic ray or fast neutral precursor. 

2. OBSERVATIONS AND RESULTS 

The spectrosopic observation was carried out on 2004 
October 1. Longslit Echelle spectroscopy of the Ha line 
was performed using order-blocking filters (HDS stan- 
dard setup "stdHa"). This gives a spectral coverage 
of 6540 A~ 6690 A over the 60" of slit length. The 
slit was centered at a(2000), 6(2000) = (00 h 25 m 56!5, 
64° 09' 28") with position- angle of 76° (measured E of 
N), covering both Knot g and its PIP simultaneously 
(Fig. Ula)). A total of 3x30 minutes of source expo- 
sure was obtained, and the same amount of exposure for 



2 



Lee et al. 



nearby sky. The spectrum was binned by 2 along the 
slit direction and 4 along the dispersion direction before 
the readout. The pixel scale after the binning is 0.27" 
pixel -1 (9.3 x 10 15 cm pixel -1 at a distance of 2.3 kpc 
(|Kamper fc van den Berehl 11971 s ) and 0.08 A respec- 
tively. The slit width was 2", which gives velocity resolu- 
tion of 17 km s -1 . The seeing was 0'.' 5 ±0'.'l. The process- 
ing of the SUBARU data included a typical CCD prepro- 
cessing (including overscan correction and flat fielding) 
and two-dimensional spectral extraction. A wavelength 
calibration solution is obtained from the spectrum of a 
Th-Ar lamp. The source spectrum was sky-subtracted, 
and normalized using the spectra of standard stars. The 
uncertainty in the wavelength calibration is estimated to 
be around 0.2 kms -1 at the wavelength of Ha. 

In Fig. 1(b), we present the fully processed two dimen- 
sional spectrum of Ha line. The bright patch at the bot- 
tom with three distinct local emission peaks corresponds 
to Knot g, and the faint emission extending to the top 
of the image corresponds to the PIP. The average Ha 
spectrum of Knot g and that of PIP arc shown together 
in Fig. [5] The velocity width of the Knot g narrow com- 
ponent line is clearly larger than that of the PIP Ha line. 
In addition, the velocity centroid of the Knot g narrow 
component is slightly redshifted (5.5 ± 0.6 kms -1 ) rela- 
tive to that of the PIP Ha line. And as clearly seen in 
Fig. H£b), the Knot g spectrum shows a very broad (~ 
1,000 kms -1 ), faint Ha line. The small peak near 900 
kms -1 is the [N II] A6583.4 A line from the PIP. 

The spectrum of the PIP is well fitted by a single Gaus- 
sian, and yields a FWHM of 33.8 ± 0.8 kms" 1 and cen- 
troid velocity of —35.8 ± 0.6 kms -1 (in LSR frame). 7 
The measured FWHM of [N II] A6583.4 A in the PIP 
is ~ 23 kms -1 . If the broadening were purely thermal, 
the widths of the lines whould be inversely proportional 
to the square root of their atomic masses. The expected 
width of the [N II] line in this case is 1/V14 = 0.267 
of Ha. As this is significantly narrower than what is 
observed, a significant amount of nonthermal broaden- 
ing is suggested. If we simply assume that the observed 
line widths are a convolution of thermal and nonther- 
mal broadenings, the estimated thermal temperature is 
- 13, 000 K, consistent with 12, 000 K of GOO. It is also 
possible that the large line width is due to a residual of 
Galactic Ha emission. 

An adequate fit to the Knot g spectrum requires three 
Gaussian components. They have velocity widths of 
45.3 ± 0.9 kms -1 (narrow), 108 ± 4 kms -1 (intermedi- 
ate) and 931 ± 55 kms -1 (broad), with central velocities 
of -30.3 ± 0.2, -25.8 ± 0.8 and 29 ± 18 kms -1 , respec- 
tively. This result confirms that the Ha narrow com- 
ponent line of Knot g is redshifted and broadened rela- 
tive to that of the PIP. The broad component (FWHM~ 
1,000 kms -1 ) should correspond to the previously re- 
ported ~ 2,000 kms -1 component (|Ghavamian et akl 
2001). The relatively narrower width of our observation 
is probably due to the insensitivity of our spectroscopic 
configuration to the very broad line, although the pos- 
sibility of temporal variation (e.g., by crossing a density 

7 throughout this paper, we give line widths corrected for instru- 
mental broadening and velocity in LSR frame. 



jump) does exist. The abrupt density increase might be 
possible if the shock is propagating int o the edge out- 
skirts of the dense cloud (|Lee et aLl l2004) . The existence 
of the intermediate width component has already been 
reported by GOO. It may be produced when slow pro- 
tons picked up by the postshock magnetic field undergo 
a secondary charge exchange. Alternatively, it might be 
an artifact of the assumption of Gaussian distributions, 
which would not necessarily be appropriate if the mo- 
tions are non-thermal. 

The characteristics of the narrow component Ha line 
of Knot g are consistent with GOO, except that the ve- 
locity centroid of the line in our data is significantly dif- 
ferent from theirs (wlsr = — 30.3± 0.2 from our data vs. 
— 53.9 ± 1.3 kms -1 from GOO). We have carefully as- 
sessed our wavelength calibration, but couldn't find any 
significant source of error. The wavelengths of night sky 
lines observed near Ha matched well with those of VLT 
UVES night-sky emission line catalog (jHanuschild l2003h 
within error of 0.02 A (~ 1 kms -1 ). The Ha spectrum of 
blank sky is also con sistent with nearby WHAM spectra 
(Haff ner et al.l [20031 ). Furthermore, the central velocity 
of [N II] A6583.4 A emission line from the PIP is con- 
sistent with that of Ha. Although we believe that our 
wavelength calibration solution is quite secure, the dis- 
crepancy with the previous observation should be con- 
firmed by an independent observation. For the rest of 
this paper, we mainly concentrate on the line width and 
the relative variation of central velocity. 

3. LOCATION OF THE SHOCK FRONT AND DISCOVERY 
OF A THIN (10 16 CM) PRECURSOR 

Fig. [3] shows the spatial variation of the Ha line profile 
along the slit. In the top panel, we separately plot the 
integrated intensities of the representative Ha of narrow 
and broad components, together with their sum (refer 
to the figure caption for the velocity range of each com- 
ponent). In the middle and bottom panels, we plot the 
central velocities and FWHMs at each slit position from 
the single Gaussian fit to the full profile. In all three 
panels, the x-axis represents the pixel offset from the 
southwest end of the slit. Although three Gaussian com- 
ponents are actually required to adequately fit the Ha 
profile of Knot g, the non-uniqueness of the multiple- 
component fitting can complicate the interpretation of 
the different emission line components. Therefore, we 
plot the fit result from a single Gaussian and describe 
our interpretation of this fit in the region of Knot g as 
following: First, the central velocities closely match those 
of the narrow component. Fitting with multiple compo- 
nents leaves the fitted line centroids unchanged to within 
the errors. Second, the FWHMs basically trace the spa- 
tial variation of the narrow component width, but the 
presence of the intermediate and the broad components 
contribute significantly to the width. When fit with three 
Gaussian components, the large widths of the Ha nar- 
row component in Knot g from the single profile models 
(60 - 80 kms -1 ) are reduced to 40 — 50 km s 

Inspection of Fig. [3] reveals that there exists only one 
location (pixel offset 27, marked as a dotted vertical line) 
where the intensity and width of the Ha line exhibits an 
abrupt jump. The jump is most noticeable for the broad 
component, where the flux is virtually zero toward the di- 
rection of ambient medium (i.e., toward positive pixel off- 
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Fig. 1. — (a) Ha image of Tycho Knot g with position of the 
2" x60" slit overlaid, (b) The fully-processed Ha 2-d spectrum 
(only the western part of the slit is shown). The bright emission 
knot at the bottom is Knot g. 

sets). The central velocity of the narrow component also 
shows rapid change around this location. The FWHM 
also steeply increases behind this point, but this might 
be largely due to the sudden appearance of the broad 
component. The fact that the intensity of the broad 
component, which is associated with the postshock gas, 
shows a significant jump at this location suggests that it 
corresponds to the location of the shock front. 

The Ha intensity in the PIP region is nearly constant 
along the slit, but shows an rapid rise just before the 
shock front, i.e., from pixel offset 32 to 27 in Fig. [3] The 
intensity of the Ha line in the PIP region is about 10% 
of the observed peak value in Knot g, and it reaches 
about half the peak very near the shock. The line width 
also increases from 30 to 45 kms -1 within this region. 
We propose that the steep increase of the Ha intensity 
accompanied with line broadening corresponds to a thin 
precursor where gas is heated and accelerated. Unlike 
the line broadening, the observed velocity centroids shift 
slightly behind the shock instead of the precursor region. 
This does not necessarily indicate that the broadening 
and the Doppler shift takes at different region as this 
is likely due to geometrical projection effects. The Ha 
intensity profile gives an e- folding thickness of 1.4 ± 0.4 
pixel (after accounting for seeing), which corresponds to 
(1.4 ± 0.4) x 10 16 cm at a distance of 2.3 kpc. 

The observed emission of Knot g shows a few local 
peaks indicating possible substructure, e.g., a collection 
of shock tangencies seen in projection along the line of 
sight. This leads to a possibility that the 'precursor' 
is simply the results of geometric projection of fainter 
Balmer-dominated filament lying just ahead of Knot g. 
If we assume that this filament has a line profile simi- 
lar to that of Knot g, i.e., have a same broad-to- narrow 
ratio, then a detectable amount of broad component is 
expected. The flux profile of broad component plotted 
in Fig. [3] would show similar gradual increase in the pre- 
cursor region, which is not seen. Also, no hint of a broad 
component is found in a summed spectra of precursor 
region, which is expected to show up. Examining the 
archival Chandra observation of Tycho's SNR (the anal- 
ysis is presented by Warren et al. 2005), we did not 
find evidence of X-ray emission extending ahead of Knot 
g. Therefore, there is no strong supporting evidence in 
favor of a projection, although we cannot rule out the 
possibility. 

4. NATURE OF THE THIN PRECURSOR 

As the shock is nearly tangential to our line of sight, 
the actual amount of bulk gas acceleration could be much 
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Fig. 2. — (a) Ha spectrum of Knot g and the PIP. The spectrum 
of Knot g is fitted with three Gaussian components, and two nar- 
rowest components are displayed (thin dashed) together with their 
sum (thin solid), (b) Same as (a) except the full observed velocity 
range is presented and the spectra are binned. The Knot g broad 
component from the above fit is displayed as a thin solid line. 

larger than the observed Doppler shift. The shock an- 
gle can be inferred from the radial velocity shift of the 
broad c omponent of Knot g r elative to the narrow com- 
ponent (jChevalier et al.lfl980f ). which is measured to be 
~ 60 kms -1 . As our observation could be insensitive to 
this broad component, usin g this value is rather con ser- 
vative. On the other hand, Ghavam ian et all (|2001l ) re- 
ported a redshift of ~ 130 kms" , but their field of view 
is different from ours. We take these two values as limits 
and give shock angle of 2° — 5° assuming a shock velocity 
of 2,000 kms - . We consider 5 kms -1 to be the repre- 
sentative redshift of the narrow component compared to 
unperturbed medium, which gives actual acceleration of 
60 - 130 kms -1 . 

The line width of the narrow component at the shock 
front (~ 45 kms -1 ) corresponds to a neutral hydrogen 
temperature of 40, 000 K, if the line broadening is purely 
thermal, or lower, if there is a significant non-thermal 
broadening such as a wave motion. Neutral hydrogen 
atoms and protons may have similar velocity distribu- 
tions due to their charge exchange interactions. But that 
of electrons can be different, which would depend on the 
heating mechanism in the precursor. In the following, 
we estimate the electron temperature in the thin precur- 
sor (TP) from the observed intensity increase of factor 
5. Since the observed spectrum is an integrated emis- 
sion along the line of sight where a significant contribu- 
tion from PIP region is expected, the actual emissivity 
increase within the precursor should be much greater. 
We assume that regions of the PIP and the TP are 
represented by two concentric shells with thicknesses of 
<5pip ~ 10 16 cm and <5tp — lpc (GOO), respectively, and 
that both have an inner diameter of the Tycho (5.4 pc). 
Then the ratio of path length through each shell along 
the tangential direction of the shock front is ~ 18. This 
implies that the emissivity increase in the TP could be 
as large as a factor of 90. This value should be regarded 



4 



Lee et al. 




-20 
I" 25 



-30 



u 

-35 

80 

„ 70 
E 60 

2 50 
I 

£ 40 



+r+, 



¥+++ ++++ 



30 



20 





, + + + V + 






+ + 






+ 














+ ^+ + + +^ + 
+ ++ ++-_+_- 

+ 



10 



20 30 40 
pixel offset 



50 



60 



Fig. 3. — (top) Integrated fluxes along the slit position of 
Ha spectrum for a given velocity range. The thin solid line 
shows emission in the range —60 < ulsr < 10 kms -1 , rep- 
resentative of Ha narrow component, while the thick solid line 
shows emission in the range —400 < ulsr < —120 kms -1 and 
+100 < i>lsr < +350 kms -1 , for the broad component. The 
thick dashed line shows the summed intensity for both components. 
(—400 < ^lsr < +350 kms -1 ). Middle and bottom panels: cen- 
tral velocities and FWHMs from fits with a single Gaussian. The 
proposed location of shock front is marked as a vertical dotted line. 

as an upper limit as it is likely that the filament is nearly 
flat and tangent to the line of sight. The collisional ex- 
citation rate of Ha at T c ~ 40, 000 K is 2,000 times 
greater than that of 12, 000 K which is the temperature 
of the PIP region (GOO). This value greatly exceeds the 
emissivity increase of 90, and implies a few possibilities. 
The electron temperature can be less than 40, 000 K, ei- 
ther if observed Ha line width has significant nonthermal 
broadening, or if T c is intrinsically less than T p . On the 



other hand, the estimated emissivity increase might be 
explained if the emissivity increase due to high T c is sup- 
pressed by ionization of neutrals. Since the existence of 
Blamer-dominated filaments requires a significant frac- 
tion of neutral hydrogen atoms to survive within precur- 
sor, this possibility is less favored. 

The two likely candidates for this precursor are fast 
neutral and CR precursors. The momentum and en- 
ergy carried by upstreaming fast neutrals can be large 
enough to explain the observed heating and accelera- 
tion in the precursor, but the available model calcula- 
tions do not predict significant hea ting by these neutrals 
(|Lim fc Ragalll99d iKorrecld 12005( 1. Furthermore, It is 
difficult to reproduce the small range of FWHMs ob- 
served for narrow component Ha lines from shocks of dif - 
ferent velocities ([Smith et al.l Il994t iHester et all 1 1 994( 1. 
The observed characteristics of the precursor are consis- 
tent with a CR precursor. C R acceleration in the shock 
does require a precursor (e.g.. lBlandford fc EichleJfl987l ) 
and there has been growing evidence of C R acceleration 
in yo ung SNRs including Tycho itself (e.g.. lWarren et aTl 
2005). Significant heating and acce leration are expected 
to ha ppen in the CR precursor fBlandfor d fc Eichled 
1987). Interaction of CR particles in the upstream gen- 
erates Alfven waves and signific ant amplification of mag- 
netic field has been suggested (|Bell fc Lu cek 2001( 1, Al- 
though a measurement of the preshock magnetic field 
is hardly available, a high value (4 /xG) of preshoc k 
magnetic field is sugg ested for Tycho (|Volk et al.112001 ). 
which may explain the line width of the narrow compo- 
nent. 

To conclude, our SUBARU observation clearly has 
shown the line broadening and the Doppler shift between 
the preshock gas and postshock gas. This strongly sug- 
gests the existence of the thin precursor. Furthermore, 
the precursor itself is likely to be resolved. Given the 
lack of observational constraints on CR ion acceleration 
in SNR shocks, the tentatively measured width (~ 10 16 
cm) of the thin precursor whose primary candidate is the 
CR precursor is promising. Follow-up observation, such 
as a high resolution imaging with HST, would clearly 
resolve the structures of the precursor. 
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